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ABSTRACT 

We report the discovery of excess 4.5 and 8/im emission from three quiescent black hole low-mass X- 
ray binaries, A 0620-00, GS 2023+338, and XTE J1118+480, and the lack of similar excess emission 
from Cen X-4. The mid-infrared emission from GS 2023+338 probably originates in the accretion 
disk. However, the excess emission from A 0620—00 and XTE J1118+480 is brighter and peaks at 
longer wavelengths, and most likely originates from circumbinary dust that is heated by the light 
of the secondary star. For these two sources, we find that the inner edges of the dust distributions 
lie near 1.7 times the binary separations, which are the minimum radie at which circumbinary disks 
would be stable against tidal disruption. The excesses are weak at 24 ^m, which implies that the 
dust does not extend beyond about 3 times the binary separations. The total masses of circumbinary 
material are between 10 22 and 10 24 g. The material could be the remains of fall-back disks produced in 
supernovae, or material from the companions injected into circumbinary orbits during mass transfer. 
Subject headings: circumstellar matter — stars: individual (A 0620—00, GS 2023+338, XTE 
J1118+480, Cen X-4) 



1. INTRODUCTION 

One of the main limiting factors in determining the 
masses of the compact objects in low-mass X-ray binaries 
is uncertainty in the fraction of the infrared light that is 
produced by the companion star (e.g., Hynes, Robinson, 
& Bitner 2005). When LMXBs are accreting at low rates, 
the optical and infrared light from these systems is dom- 
inated by the companion stars, so their mass functions 
can be measured to high accur acy based on the Doppler 
motio n of the companion (e.s.. Mc Clintock fc Rem illard 
2006) . The inclination then must be constrained by mod- 
eling the modulations in the optical and infrared light 
curves that are produced by the varying aspect of the 
distorted, Roche-lobe-filling companion. Generally, the 
contribution of the accretion disk to the optical and near- 
infrared emission is uncertain, and varying the fraction 
of the light that is assumed to be produced by the ac- 
cretion disk can lead to differences of a factor of 2 in the 
derived mass of a compact object (e.g., Gclino, Harrison, 
& Orosz 2001). Fortunately, the broad-band spectrum of 
a multi-temperature accretion disk is significantly flatter 
than that of a stellar photosphere, so mid-infrared obser- 
vations could constrain the relative contributions of the 
two components. 

However, the spectra of LMXBs in the mid-infrared 
have not been well studied, and several indirect lines of 
evidence suggest that these systems might contain cir- 
cumbinary material that could emit in the mid-infrared. 
gFirst, some of the white dwarf analogs to LMXBs, 
cataclysmic variables (CVs), exhibit spectral features 
that lie at the mean radial velocities of the systems 
( Solh emTfcSior] ITiwl and excess mid-infrared emis- 
sion l|Howell et alJ 120061) that could be in terpreted as 
n circumbinary material (see also ITaam et al.l 
Dub us et all 12004^ . Second, the supernovae that 
produced the compact objects could have left fall-back 
disks around the binaries. Indeed, the first fall-back 
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disk has recently been found around a young, highly- 
magnetized neutron star (Wang, Chakrabarty, & Ka- 
plan 2006). Third, the planets around the isolated mil- 
lisecond pulsar PSR 1257+12 ijWolszczan fc Frailll99^ 
could not have survived the supernova that produced the 
neutron star, and must h ave formed afterward (see, e.g., 
IMiller fc Hamiltonl2001|) . Millisecond pulsars are usually 
assumed to have been spun up by accretion as LMXBs, in 
which case planets could form from material present dur- 
ing the binary phase. Therefore, to search for evidence 
of circumbinary material, we have observed four nearby, 
quiescent LMXBs with the Spitzer Space Telescope. 



2. OBSERVATIONS AND DATA ANALYSIS 

We chose the LMXBs in our sample to be detectable 
with Spitzer if they contained optically-thick circumbi- 
nary _jlj^sjDa^ivel^jruHnm the mass donor stars 
(see lJurall2003l ITaam et al.ll2003l and below). Based on 
a simple model, we chose sources with: (1) K magni- 
tudes brighter than 17, (2) locations more than 50° in 
projection from the Galactic center, and (3) no 2MASS 
sources within 5" that were brighter than our targets. In 
Tabled we list the positions, the orbital periods, esti- 
mates of the primary masses, the spectral types of the 
companions, and the quiescent K magnitudes of the four 
systems in our sample (see the table notes for references) . 

Our measurements were taken in the 4.5 and 8.0 /im 
bands with Infrared Array Camera (IRAC), and in the 
24 /jm band with the Multiband Imaging Photometer 
(MIPS; Tab. We used the post-basic-calibration data 
provided by the Spitzer Science Center (SSC) for most 
of our analysis. However, the MIPS image of A 0620—00 
contained latent features with a low spatial frequency 
and a ~2% amplitude that were left by a previous ob- 
servation of a bright, extended source. We corrected the 
image by creating a flat field from the median of the 
individual dithered images, dividing each snapshot by 
the flat, and re-creating the mosaicked image using the 
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Fig. 1. — Mid-infrared images of the LMXBs in this survey. The 
red channel is 24 /im, the green 8.0 fim, and the blue 4.5 (im. The 
circles denote the positions of the optical counterparts. 



script provided by the SSC. 2 Three-color images centered 
on each LMXB are displayed in Figure ^ Each target 
is detected at 4.5 and 8.0 /im. Only A 0620-00 and 
GS 2023+338 are also detected at 24 /im. We com- 
puted the IRAC fluxes of each source using the point- 
spread-function-fitting routine APEX from the SSC, and 
the MIPS fluxes and upper limits using aperature pho- 
tometry (Tab. El and Fig. EJ) . 

To understand the origin of the mid-infrared emis- 
sion, in Figure [21 we plot for each LMXB the observed 
and de-reddened fluxes in the infrared and optical bands 
(from Tab. El and the references in Tab. The opti- 
cal and near-infrared intensities of quiescent LMXBs of- 
ten vary by several tenths of a ma gnitude on time scales 
of years (not coun ting outbursts: iChevalier et ail 119891: 
iGelino et al.l |2001). presumably because of changes in 
the accretion flow. None of the fluxes were obtained 
simultaneously, so when comparing the flux measure- 
ments in Figure [21 one can only expect them to be self- 
consistent to within «30%. To estimate the contribu- 
tions of the companions to the spectra, we have ob- 
tained model stellar spectra computed with the next- 
generation Phoenix code (kindly provid ed by T. Barman; 
lHauschildt. Lowenthal. fc Baronl 120011) that correspond 
to the temperatures and radii of the spectral types in 
Tabled We have estimated the contribution of the ac- 
cretion disk using a standard Sh akura fc Sunvaevl (jl973f) 
model, assuming the disk extends out to the tidal trunca- 
tion radius (<60% of the binary separation; Frank, King, 
& Raine 1992), is inclined by 60° to our line of sight, and 
has the largest accretion rate that is consistent with the 
optical photometry (see the caption of Fig. [21 for val- 
ues) . We have confirmed that our conclusions are robust 
against changing the assumed spectral types of the com- 
panions within the values reported in the literature, and 
against varying the parameters of the accretion disk so 
long as the models are consistent with the optical pho- 
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Fig. 2. — The spectra of the LMXBs in our sample. The filled 
circles are the observed fluxes in each band, and the open circles 
are the de-reddened fluxes (see notes to Tab. 1 for references). 
The upper limits on the 24/mi fluxes for XTE J1118+480 and Cen 
X-4 are indicated with arrows. The uncertainties in the observed 
fluxes are also indicated, but are generally smaller than the symbol 
sizes. The solid grey lines are synthetic stellar spectra for the 
spectral type in Table 1. The dotted grey lines are model spectra 
of viscously-heated accretion disks, with accretion rates of 3 X 10 14 
g s" 1 for A 0620-00, 7 x 10 16 g s" 1 for GS 2023+338, 4 x 10 13 
g s" 1 for XTE J1118+480, and 4 X 10 14 g s" 1 for Cen X-4. The 
grey dashed lines illustrate the spectra of blackbodies used in the 
text to estimate the solid angles of the material responsible for the 
mid-infrared excesses from A 0620-00 and XTE J1118+480. The 
black lines are the sum of these model components, normalized to 
match the de-reddened K flux. 



tometry. 

We find that the spectrum of Cen X-4 is consistent 
with that expected for the companion star, with a pos- 
sible contribution at long wavelengths from the accre- 
tion disk. There are clear excesses of mid-infrared flux 
above that expected from the companions of A 0620- 
00, GS 2023+338, and XTE Jl 118+480. In order to 
conservatively estimate the amounts of the mid-infrared 
fluxes that are not produced by the companions, we com- 
puted the maximum contributions of the companions' 
photospheres by normalizing the model stellar spectra to 
match the de-reddened K fluxes, and reported the dif- 
ference between the predicted and observed mid-infrared 
fluxes as the excesses in Tabled (AS\). 

3. DISCUSSION 

Several LMXBs have now been detected in outburst 
in the mid-infrared (e.g . . Smith, Beall, fc Swain 1990 ; 
I van Paradiis et alJll990t lHarrison fc Oehrdll99ll Il994t 
Migliari et al. 2006), but these are the first detections of 
quiescent systems. In outburst, the mid-infrared emis- 
sion could originate from an accretion disk in which 
M > 10 18 g s _1 l|Smith et al.lll990|) . from free-free emis- 
sion in a strong wind d riven from the accretion disk 
l)van Paradiis et aDll99f)D . or from radio jets with flat 
spectra extending t o infrared wav elengt hs l)Fenderll200ft 
iHoman et al] 120051 iMigliari et alJ 120061) . We use these 
hypotheses as a starting point for understanding the mid- 
infrared emission from our sample of quiescent LMXBs. 

The mid-infrared excess from GS 2023+338 could orig- 
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TABLE 1 
Target Quiescent LMXBs 



Source 


Optical 


RA DEC 


D 


-^"orb 


Mi 


Sp. Type 


A v 


K 






(J2000, degrees) 


(kpc) 


(hr) 


Mq 








A 0620-00 


V616 Mon 


95.68561 -0.345628 


1.2+0.1 


7.8 


8.7-12.9 


K4V 


1.2 


14.55(6) 


GS 2023+338 


V404 Cyg 


306.01594 +33.86728 


2.2-3.7 


155.3 


10.1-13.4 


K0III 


4.0 


12.50(5) 


XTE J1118+480 


KV UMa 


169.54498 +48.03678 


1.8+0.5 


4.1 


6.5-7.2 


K6V 


0.06 


16.9(2) 


Cen X-4 


V822 Cen 


224.59135 -31.66872 


1.2 


15.1 


1.4 a 


K5V 


0.3 


14.66(8) 



Note. — The positions ar e taken from 2MASS, and the distances, orbital periods, a nd mass fun c tion f rom 
IMcClintock fc Remillardl 120061) . For the remaining d ata, the references ar e — A 0620-00: IGelino et all 120011) for 
all photometry and the spectral type — GS 2023+338: ICasares et alj 12003) for all photometry and t he spectral type 
— XTE J1118+480: spe ctral type from lTorres et alj 120041): extinction fro m IMcClintock et all 120031) : J and K from 
IMikoIaiewsa et alJ 120051) ; B from USNO; an d R from lZurita et alJ 120031) — Cen X-4 : spectral type from Shahbaz, 
Naylor, & Charles ( 1993): ITbrres et alJ i2002T) : J, H, K from 2MASS; other bands from lShahbaz et all 119931) . 

a We assumed this mass for the neutron star based on those measured for binary radio pulsars ( Thorsctt & Chak rabartvl 
1999). 



TABLE 2 

Spitzer Measurements of Quiescent LMXBs 



IRAC MIPS 



Source 


Date IRAC 


Tirac 


54.5 


AS4.5 


Ss.o 


A5 8 .o 


Date MIPS 


Tmips 


S24 


A5 24 






GO 


(A«Jy) 


C#*Jy) 


(A*Jy) 


(A«Jy) 




(s) 


(i*Jy) 


(A*Jy) 


A 0620-00 


2005 Mar 25 


400 


448(13) 


194 


249(10) 


149 


2005 Mar 06 


180 


54(18) 


43 


GS 2023+338 


2004 Oct 09 


36 


3020(90) 


670 


1450(40) 


500 


2004 Oct 16 


30 


153(70) 


46 


XTE J1118+480 


2004 Nov 21 


400 


46(1) 


17 


45(7) 


34 


2005 May 13 


240 


<16 


<16 


Cen X-4 


2004 Aug 12 


300 


199(6) 


-90 


95(17) 


-14 


2005 Aug 28 


150 


<30 


<30 



Note. — Wc list the date and duration of each observation (the 4.5 and S /im data were obtained in parallel), the fluxes measured 
(S\, where A is the central wavelength of each band), and estimates of the minimum amount of excess mid-infrared emission above 
that expected from the photospheres of the companions (AS\; see text). The quoted uncertainties are the larger of the standard 
deviation of the fluxes in the multiple frames wc obtained, or the systematic uncertainty in the absolute calibration (3% for IRAC, 
and 5% for MIPS). The negative values for Cen X-4 reflect the fact that extrapolating a A 2 spectrum expected for the Raylcigh- Jeans 
tail of the photosphere from the K-band flux predicts more mid-infrared flux than is observed; this is likely because the infrared flux 
is variable. 



inate in an accretion disk, because the 2.2-24 fixa fluxes 
in Figure [21 follow the v~ 2 law that one would expect 
from the Rayleigh- Jeans tail of a black body. We find 
that our model for a vis cously-heated accretion disk 
( Shak ura fc Sunvaevlll973(l from §2 adequately describes 
the excess mid-infrared flux from GS 2023+338 in Fig- 
ure 

In contrast, the 4-24 fim fluxes from A 0620—00 and 
XTE J1118+480 do not follow the v~ 2 law expected for 
the Rayleigh-Jeans tail of a blackbody, because there is 
too much flux at 8 /im. A blackbody spectrum that 
peaks near 8 fjxa would have a temperature of only 640 
K. For these two sources, we plot blackbody spectra that 
match the peak of the mid-infrared excesses (ASg.o) m 
Figure 12 The sum of the spectra of stellar photospheres 
and optically-thick blackbodies match the observed mid- 
infrared fluxes well. However, the inferred emitting ar- 
eas are large compared to the binary separations. For a 
Planck spectrum, the solid angles of the emitting regions 
correspond to circular radii of 



R = 4 x 10 



11 / ASgpm \ 



1/2 



D 
kpc 



cm, 



(1) 



where D is the distance to the source. In contrast, the 
disks will be contained within <60% of the binary sepa- 
rations a, which are given by: 



3xlO U MV3(l +g )V3p2/3 cmj 



(2) 



where Pday is the orbital period in days, M is the mass 
of the accretor in s olar masses, and q is the mass ratio 
of the binary (e.g.. lFrank et al1l!992|) . Given the fluxes 
in Table |21 we find that the excess mid-infrared emission 
originates from regions «2 times larger than the orbital 
separations of GS 2023+338 and XTE J1118+480, and 
>4 times larger than the accretion disks. Therefore, the 
mid-infrared emission is produced in a region that ex- 
tends beyond the binary orbit. 

Jets could produce an emitting region larger than the 
binary orbit, but given the spectra of the mid-infrared 
excesses we suggest that jets are only minor contribu- 
tors. Two LMXBs in our sample have recently been de- 
tected as radio sources in quiescence: GS 2023+338 with 
flat-spectrum radio emission with an i ntensity of 350/iJ y 
over the frequency range 1.4-8.4 GHz ijGallo et all 20051. 
and A 0620-00 w ith a radio flux of 50 /xJy at 8.4 GHz 
l|Gallo et alJl2006D . The 24/im excesses from A 0620-00 
is equal to the flux in the radio, so it could be produced 
by a jet with a flat (S u cx v a , with a w 0) spectrum 
between the radio and mid-infrared. For GS 2023+338, 
the 24 /im flux is slightly lower than the radio flux, which 
would imply that its jet has a steeper a > —0.3 spec- 
trum. However, the 8 /urn excesses from A 0620—00, GS 
2023+338, and XTE J1118+480 are not consistent with 
a flat-spectrum radio jet, because for XTE J1I18+480 
and A 0262—00 they are 4 times larger than the excesses 
at 24 /im, and for GS 2023+338 it is 10 times larger. If 
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we assume the 24 /xm emission is from a flat-spectrum 
jet, then it contributes < 25% to the 8 /im flux from 
XTE J1118+480 and A 0620-00, and <10% from GS 
2023+338. 

We suggest that the excesses from A 0620—00 and 
XTE J1118+480 originate from circumbinary dust that 
re-processes the light of the companions. To estimate 
the masses of the dust, we assume it is contained in 
optically thick disks. The disks could lie as close to 
the center of masses of the binaries as 1.7a, at which 
point they would b e tidally truncated l)Taam et al]l2003t 
iDubus et alJ 12004) . We can estimate the temperature 
profile of such disks by assuming they are flat, in which 
case 

where Tdisk and i?disk are the temperatures and radii 
of the disks, and T» and i?* ar e the temperatures and 
radii of th e companions (e.g., R uden fc Pollackl 119911: 
I.Tural 12003ft . Using Equation and the parameters in 
Table ^ we would expect the inner edges of the cir- 
cumbinary disks to have temperatures of «600 K, which 
is consistent with our detection of excesses that peak 
near 8 /im. The lack of excess flux at 24/xm implies that 
material does not re-process much stellar light beyond 
radii «3 times larger than the binary separations. If we 
assume the disks are compose d of dust with an o pacity 
X„~300 cm 2 g _1 at 8 ^m fe.g.. lDraine fc Leell984D . then 
for them to have an optical depth of r~l, the disks would 
only need to contain ~10 22 g of dust. If we assume the 
emitting material is optically thin, we find dust masses 
th at are similar wi thin a factor of a few (e.g., Eq. 3 
in IE vans et al]ll997l) . If the gas-to-dust ratio is similar 
to that of the interstellar medium, ^100, then the total 
mass of the circumbinary material could be ~10 g. 

4. CONCLUSIONS 

The detection of excess mid-infrared flux from A 
0620-00 and XTE Jl 118+480 (Fig. |2J provides evi- 
dence that circumbinary material is present around some 
LMXBs, but its spectrum suggests that it will only be de- 
tectable in the mid-infrared. The excess emission peaks 
at 8 /xm, and if we model it as a single-temperature black- 
body, we predict that the circumbinary material pro- 
duced negligible flux shortward of «3 /im. Moreover, the 
lack excess mid-infrared emission from Cen X-4 demon- 
strates that such material is not ubiquitous. Finally, for 
GS 2023+338, the mid- infrared emission appears to orig- 
inate from a hot accretion disk, and the contribution of 
the accretion disk to the mid-infrared light can be pre- 
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dieted based on the optical measurements (the dotted 
line in Fig. EJ). These results provide reassurance that 
efforts to determine the inclinations of LMXBs by mod- 
eling the ellipsoidal modulations of their infrared and 
optical light curves are not compromised by the presence 
of circumbinary material. 

For A 0620-00 and XTE Jl 118+480, we suggest that 
the circumbinary disks are either the remains of fall-back 
disks produced in the supernoyae that formed the com- 
pact objects Ce.g.. lWang et al.ll2006t iCordes fe Shannonl 
2006), or material injected into circumbinary orbits dur- 
ing the proces s of mass loss by the Roche-lobe filling com- 
panions (e.s.. iDubus et aT]l2004fl . These circumbinary 
disks contain ^10 22 g of dust, which for a standard dust- 
to-mass ratio of ^100 implies a total mass of ~10 24 g. If 
the disks are produced by matter ejected into circumbi- 
nary orbits during the process of mass transfer, then it 
represents only a tiny fraction of the ~10 32 g (~0.1 Mq) 
that will be lost by the comp anions over the lifet i mes o f 
these LMXBs. Alternatively, iGordes fc Shannonl (2006) 
have proposed that fall-back disks containing >10 24 g 
of material could form asteroids and inject them into 
the magnetospheres of pulsars with a high enough rate 
to explain the the observed intermittency in their ra- 
dio pulse s. Asteroids might not su rvive the outbursts of 
LMXBs l)Miller fc Hamiltodl27)0l|) . but the mid- infrared 
excesses that we have identified could be the remnants 
of similar fall-back disks. 

However, the circumbinary material is not massive 
enough either to affect the evolution of the orbital angu- 
lar m omentum of LMXBs (e.g., ~10 29 g in T aam et alJ 
20031) . or to form pla netesimals (e.g., ~10 28 g in 
Miller & Ham ilton! 1200 lfl . As the formation and evolu- 
tion of circumbinary matter around LMXBs and CVs are 
considered further, it may be that the paucity of matter 
around A 0620-00 and XTE Jl 118+480 (and the appar- 
ent lack of circumbinary matter around the one neutron 
star LMXB in our sample, Cen X-4) may be the most 
lasting aspect of this result. 
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